Theory of the 1-point PDF for incompressible Navier-Stokes fluids 
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Fundamental aspects of fluid dynamics are related to construction of statistical models for in- 
compressible Navier-Stokes fluids. The latter can be considered either deterministic or stochastic, 
respectively for regular or turbulent flows. In this work we claim that a possible statistical for- 
mulation of this type can be achieved by means of the 1-point (local) velocity-space probability 
density function (PDF, /i) to be determined in the framework of the so-called inverse kinetic theory 
(IKT). There are several important consequences of the theory. These include, in particular, the 
characterization of the initial PDF [for the statistical model {/i, T}] . This is found to be generally 
non-Gaussian PDF, even in the case of flows which are regular at the initial time. Moreover, both 
for regular and turbulent flows, its time evolution is provided by a Liouville equation, while the 
corresponding Liouville operator is found to depend only on a finite number of velocity moments 
of the same PDF. Hence, its time evolution depends (functionally) solely on the same PDF. In 
addition, the statistical model here developed determines uniquely both the initial condition and 
the time evolution of f\. As a basic implication, the theory allows the exact construction of the 
corresponding statistical equation for the stochastic-averaged PDF and the unique representation of 
the multi-point PDF's solely in terms of the 1-point PDF. As an example, the case of the reduced 
2-point PDF's, usually adopted for the statistical description of NS turbulence, is considered. 

PACS numbers: 47.27.Ak, 47.27.eb, 47.27.ed 



1 - INTRODUCTION 

The description of fluids, and more generally of con- 
tinua, is based on the introduction of a suitable set of 
fluid fields {Z} which define the state of each fluid and 
obey, by assumption, a well-posed set of PDE's denoted 
as fluid equations. The formulation appropriate for an in- 
compressible Navier-Stokes (NS) fluid - based on the so- 
called incompressible NS equations (INSE) - is recalled 
in Appendix A, together with the basic notations and 
definitions here adopted. The fluid equations, the fluid 
fields and the related initial and boundary conditions, 
are considered either as deterministic or stochastic, re- 
spectively for regular or turbulent flows (see Appendix 
A, Subsections A. 2 and A. 3). 

The statistical treatment of fluids usually adopted for 
turbulent flows (which may be invoked, however, to de- 
scribe also regular flows) consists, instead, in the intro- 
duction of appropriate axiomatic approaches denoted as 
statistical models. These are sets {/, T} formed by a 
suitable probability density function (PDF) and a phase- 
space T (subset of W 1 ) on which / is defined. By defini- 
tion, a statistical model {/, T} of this type must permit 
the representation, via a suitable mapping 



{f,T}^{Z}. 



(1) 



of the complete set (or more generally only of a subset) 
of the fluid fields {Z} = {Zi, i = 1, n} which define the 
state of the same fluids. Their construction involves, be- 
sides the specification of the phase space (r) and the 
probability density function (PDF) /, the identification 
of the functional class to which / must belong, denoted 



as {/} . As a consequence, the fluid fields Z% £ {Z} are 
expressed in terms of suitable functionals (called mo- 
ments) of /. In the case of an incompressible (and isen- 
tropic) Navier-Stokes fluid (INSF) the latter are iden- 
tified with {Z} = {V,p, St}, V and p indicating re- 
spectively the fluid velocity, the fluid pressure, both as- 
sumed strong solutions the incompressible Navier-Stokes 
equations (INSE), and St the constant thermodynamic 
entropy. As an alternative, p can be replaced by the ki- 
netic pressure p\ , to be suitably defined [see subsequent 
Eq.®]. 

Goal of this paper is to propose a statistical model for 
incompressible NS fluids, described by INSE, which is 
simultaneously (a) unique, (b) statistically complete and 
(c) closed. While the precise meaning of these statements, 
in particular the relationship fT]). will be discussed in 
detail below, we anticipate that: 

• the first condition (uniqueness) imposes that the 
same PDF and its time evolution should both be 
uniquely prescribed in terms of conditional observ- 
ables (see related definitions in the Appendix A); 

• the second one (statistical completeness; see subse- 
quent Subsection ID) involves the assumption that 
the PDF should determine uniquely also a pre- 
scribed set of physical observable (or rather con- 
ditional observables), to be suitably defined. The 
latter include in particular the complete set of fluid 
fields describing the state of the fluid; 

• the third one (closure) includes both moment- and 
kinetic- closure conditions (the corresponding defi- 
nitions are given again in Subsection ID), involving 
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- respectively - the requirements that there exists a 
closed set of fluid equations and that the statistical 
equation advancing in time the PDF depends only 
the same PDF, via a finite number of moments (of 
the PDF). 



1A - CSM-inspired models for the 1— point Liouville 
equation 

A well-known example of statistical model for incom- 
pressible NS fluids is provided by the so-called statisti- 
cal hydromechanics developed originally by Hopf l| and 
later extended by Rosen [2[ and Edwards [3( (HRE ap- 
proach). This relies on the introduction of the 1— point 
(or local) velocity-space PDF, /i, to be intended as the 
conditional PDF of the velocity v (kinetic velocity) with 
respect the remaining variables. In the HRE approach 
these are identified with (r,i), where (r,t) 6 f2 X I, while 
fi = /i(r, u,t;Z), with u = v — V(r,<) the relative ki- 
netic velocity, is identified with 



f H = 5(v-V(r,t)) 



(2) 



(deterministic PDF), fu denoting the three-dimensional 
Dirac delta defined in the velocity space U = R 3 , with v 
belonging to U = M. 3 and (r,t) G x / (with denoting 
the closure of the configuration domain fl C M 3 ). Hence, 
Jh is defined by assumption on the set spanned by the 
state vector x =(r,v), i.e., T = Q, X U, with T denoting 
the closure of the restricted phase space 



T = nxU. 



(3) 



R follows that, in this case, only the first two (velocity) 
moments of /i, corresponding to G = l,v, are actually 
prescribed in terms of the fluid fields and read 



d 3 vGf H (r,u,t;Z) = l,V(r,t). 



(4) 



Goal of the HRE approach in the case of turbulent flows 
is actually to predict 



(h(r,VL,t;Z)) = (f H (r,u,t;Z)) 



(5) 



the brackets " (•} " denoting an ensemble average, to be 
suitably prescribed, over the possible realizations of the 
fluid [a]. Rs definition, however, is not unique. For exam- 
ple, in the case of the so-called stationary and isotropic 
turbulence [3| this is usually defined so that - by assump- 
tion - it commutes with all the differential and integral 

operators (respectively, Jj,V, V 2 and J d 3 r, J d 3 v and 

n u 

J d 6 x) appearing in the NS operator [see Eq. (|94j) . in Ap- 
r 

pendix A]. As an alternative, as discussed elsewhere [HI], 
(•} can also be identified with the stochastic- averaging 
operator (|106|) defined in Appendix A. 



Nevertheless, in the HRE approach (Jh) is not di- 
rectly determined. Rather, it is replaced by a suitable 
functional of (fn) > 4>i which obeys a suitably-prescribed 
functional-differential equation (the so-called "<f> equa- 
tion" [1]). The problem of the construction of an equiv- 
alent evolution (or so-called "transport") equation for 
(Jh) has been investigated by several authors (see, for 
example, Dopazo [f| and Pope jy]). The construction of 
its formal solution is due to Monin Q and Lundgren [§] 
(ML approach; see also Monin and Yaglom, 1975 [9j and 
therein cited references). However, alternative (approxi- 
mate) statistical approaches are known, such as the GLM 
(generalized Langevin model; due to Pope [HI]). All of 
them typically rely on the assumption of the existence 
of a suitable underlying phase-space classical dynamical 
system which evolves in time the state vector x = (r, v) , 
namely the flow 

T to ,t ■ x Q -> x(t) = T to , t x (6) 

generated by an initial value problem of the type 

IE = X(x,t;Z 



X(t ) 



(7) 



Here the notation is standard [HI, [12]. Thus, T to ,t is the 
evolution operator generated by 



X(x,t;Z) = {v,F}, 



(8) 



x(x D , t D , t) being the solution of the initial- value problem 
and F(x, t; Z) a vector field (denoted as mean-field force) 
to be suitably prescribed. The definitions of the PDF 
and of its functional class depend on the type of 

relationship established between the fluid fields and the 
PDF, to be prescribed in some suitable sense. In CSM- 
inspired statistical models this is realized by means of 
a PDF, fi, which is assumed to satisfy the correspond- 
ing 1— point Liouville equation - here denoted as inverse 
kinetic equation (IKE [111 Il2|) - of the form 



L(r,v,t;/ 1 )/ 1 (r,u,i;/ 1 ) = 0. 



(9) 



Here L(r, v, t; /i) denotes the Liouville streaming opera- 
tor 

L(r, v, t; f{). = -•+—• {X(x, t; Z)-} = (10) 
d d d r „. „ , , 

while the vector field F is generally to be assumed func- 
tionally dependent on f\. In particular, in the HRE and 
ML approaches this is obtained by invoking the position 



1 H- 



ill) 



with Fh denoting the total fluid force density acting on 
each fluid element [see Eq. ([96]) in the Appendix] and iden- 
tifying fi with the particular solution 



/i — Ih 



(12) 
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[with fjj defined by Eq.©]. Hence in this case {fi} 
is manifestly the functional class of distributions of the 
form (J2J. 

The HRE and ML approaches, both characterized by 
the same statistical model {/#, T} (HRE-ML statistical 
model), belong actually to a more general class of statis- 
tical models inspired by Classical Statistical Mechanics 
(CSM) (see also Sec.2). 

In particular, this concerns the representation of all the 
relevant dynamical variables in terms of hidden variables 
Id | (see Appendix A). By definition they denote a 
suitable set of independent variables a = {o^, i = 1, k} S 
V a C R fc , with k > 1, which cannot be known deter- 
ministically, i.e., are not observable. In the context of 
turbulence theory these variables are necessarily stochas- 
tic. This means that they are characterized by a suitable 
stochastic probability density g defined on V a (sec defini- 
tions and related discussion in Appendix A, Subsection 
A. 2), while the ensemble average (•) [defined in Eq.©] 
can be identified with the stochastic- averaging (•) de- 
fined by Eq. (|106p [see Appendix A]. Hence, for turbu- 
lent flows the fluid fields - together with the PDF f\ and 
the vector field F(x, t; Z) appearing in Eqs.((8]) and ([9]) 
admit a representation of the form 



m Kqs . 

HQ 



{Z} = {Z(r,t,a)} 
fi = A(r,u,i,a;Z) 
F = F(x,t,a;Z) 



(13) 



to be defined in terms of a set of hidden variables a and a 
stochastic model {g, V a } (see again Appendix A, Subsec- 
tion A. 2). Hence, {Z} , f\ and F(x, t; Z) are necessarily 
non-observable. Nevertheless, if we assume that the fluid 
fields {Z} are uniquely-prescribed ordinary functions of 
(x,t, a) defined for all (x,t, a) 6 T x / x V a , it follows that 
they can still be considered conditional observables (see 
Appendix A, Subsection A.l). Similar conclusions apply 
to /i, and to the vector field F(x, t; Z) as well. 



IB - The Closure problem 

Based on these requirements [the positions defined by 
Eqs. fTTj) and (|12p j. in turbulence theory, i.e., when the 
fluid fields {Z} are considered as stochastic functions, 
usually the statistical description involves the construc- 
tion of an infinite set of continuous many-point PDF's 
which obey a hierarchy of statistical equations, the so- 
called ML (Monin-Lundgren 0, HI) hierarchy. There- 
fore in this case the statistical model actually involves 
the identification of {/} with the functional class of the 
many-point PDF's. As proven by Hosokawa [HI], the 
problem can be formulated in an equivalent way in the 
framework of the HRE approach, yielding the well-known 
Hopf <j) functional-differential equation. 



To date, the search of possible exact "closure condi- 
tions" for the ML hierarchy (or closure problem) - or 
equivalent, of exact solutions of the HRE approach - re- 
mains one of the major unsolved theoretical problems 
in fluid dynamics. For the ML-approach, this involves 
in principle the construction of statistical models which 
should be characterized by a finite number of (multi- 
point) PDF's, i.e., determined in such a way that the 
time evolution of the fluid fields can be uniquely deter- 
mined in terms of them. In practice the program of con- 
structing theories of this type, and holding for arbitrary 
fluid fields, is still open due to the difficulty of preserv- 
ing the full consistency with the fluid equations. In fact, 
it is well known that many of the customary statistical 
models adopted in turbulence theory - which are based 
on closure conditions of various type (see for example 
Monin and Yaglom [§] 1975 and Pope, 2000 - typi- 
cally reproduce at most only in some approximate (i.e., 
asymptotic) sense the fluid equations. 

In particular, an interesting related issue is that posed 
by the determination of the form of the 1— point PDF. 
In decaying isotropic turbulence, according to some au- 
thors (see in particular Batchelor 0]) this is predicted as 
almost-Gaussian. Although others, based on the adop- 
tion of suitable model dynamical systems for NS turbu- 
lence (Falkovich and Lebedev [l6| and Li and Menevau 
[l3]), have pointed that the tails of the 1-point PDF 
might exhibit a strongly non- Gaussian behavior, accord- 
ing to more recent investigations (Hosokawa [18[) there 
seems to be still insufficient experimental evidence for a 
generalized behavior of this type, at least in the case of 
homogeneous turbulence. 

The basic difficulty is, however, related to the proper 
formulation of a rigorous theory for the 1-point PDF 
holding for arbitrary NS fluids. In this reference, in par- 
ticular, an important issue is related to the quest of a pos- 
sible exact statistical evolution equation for the 1— point 
PDF [l4j |. holding both for deterministic and stochastic 
fluid equations, which is capable of yielding the correct 
fluid fields and holds for arbitrary initial and boundary 
conditions of the relevant related physical observables 
(or, respectively, conditional observables). 

This refers, in particular, to the subset of statistical 
models {/, T} in which the PDF / is considered as an 
ordinary function. 

An example is provided by IKT 11, 12| in which the 
1— point PDF obeys by construction, both for regular 
fluids (i.e., deterministic) and turbulent (i.e., stochastic) 
fluid fields, an IKE of the form $9§. 

This result raises the interesting question whether, in 
some suitable setting, i.e., for appropriate statistical mod- 
els and in particular in the case in which {/} is a set of 
ordinary functions, the closure problem can actually be 
solved. In this regard, there are actually two possible 
routes which seem currently available: 
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• the first route: for a prescribed statistical model 
{/, T} it involves the search of possible finite subset 
of statistical equations (formed by the s equations 
for the PDF's, /„, having n < s and s finite) which 
define a closed set of equations; 

• the second route: is based on the search of possi- 
ble alternative statistical models {/, T} , based on 
the hidden-variable representation (| 1 3 j) and math- 
ematically equivalent to the complete set of fluid 
equations for INSF, for which T coincides with the 
restricted phase-space and the PDF, identified with 
the 1— point (or local) PDF fx, satisfies an IKE 
of the type defined by (0|) 7 in which the mean-field 
force F depends functionally only from fx , via suit- 
able moments of the PDF. Statistical models of this 
type (in which the evolution operator depends only 
on fx, via its moments) are usually said to satisfy 
a kinetic closure condition. 



An important preliminary task to accomplish is to es- 
tablish the relationship of fx with the fluid fields. More 
precisely, here we wish: 

A) to assess whether - besides the complete set of fluid 
fields - the PDF is possibly related to additional observ- 
ables or conditional observables. This information, in 
fact, might provide an effective constraint on the func- 
tional class of the initial PDF - evaluated at the initial 
time (t = t a ) - {fi} ro ; 

B) to determine whether, for suitable fluid fields, there 
exists a deterministic limit for fx, namely for which it 
is a Dirac delta of the type \12) . Such a limit would 
manifestly provide a constraint on the functional class 
{/i}- 

Let us point out that for both problems a simple solu- 
tion actually exists. 



lC.a - The 1— point velocity-frequency density function 



The first approach, by far the most popular one in the 
literature and adopted by several authors for the con- 
struction of approximate closure models of the ML hier- 
archy, poses nevertheless - as indicated above - a problem 
of formidable difficulty. 

Instead, a possible candidate for statistical models of 
the second type, adopting the hidden-variable represen- 
tation (Ti3)) . is already known [LJ, [lj|. It is provided by 
the 1— point PDF fx which characterizes the IKT for an 
incompressible NS fluid 11 , 12] ■ In this case the PDF can 



always be required to be also a velocity-space probability 

density, i.e. to satisfy the normalization J d 3 v fx = 1. 

u 

In both cases, however, it must be stressed that 'a 
priori' the definition of relevant statistical model still re- 
mains essentially arbitrary. This concerns, besides the 
choice of the phase-space T (only for the first route) and 
of the functional class {fx} to which the PDF belongs, 
also the definition of the set of moments to be associated 
to the fluid fields {Z}. In particular, in IKT in principle 
the definition of fx is non-unique because higher-order 
moments of fx may be still undetermined. In addition 
fx - even if assumed as an ordinary function - still remains 
by definition non-observable. 



1C - Physical realizability conditions on the 1— point 

PDF 



In reference to statistical models based on /i,such as 
the IKT model, a natural question arises, i.e., whether 
the arbitrariness in their definition can be used, by proper 
prescription on its functional class {fx}, to determine 
it uniquely consistent not only with INSE but also with 
the relevant physical observables (or conditional observ- 
ables). In such a case fx might, in particular, be viewed 
as a conditional observable too. 



In particular, regarding the observable, here we claim 
that it can be identified with the configuration- space av- 
erage of 1— point PDF, fx(t) = /i(r, v,i, a; Z), namely 
(fx(t))n , (-)q denoting the f2— averaging operator. For a 
generic phase function a(r, v,t, a) its configuration-space 
average on SI, can be identified respectively either with 
the continuous or discrete operators 



< a(r ' V '*' a)) « = ^ 



d 3 ra(r,v,t, a) (14) 



a(r,v,t,a)} n = lirn^^ ^ a(r,v,t,a). (15) 



i=l, N 



In practice, for actual comparisons with experimental 
data, both operators can be conveniently replaced by a 
finite summation of the form 

(a(r,v,t,a)) n ^ a (A) (v,t,a) = ^ a(r. ( , v,t, a), 

* i=l,N, 

(16) 

denoting a suitable integer to be considered ^ 
1. Thus, (/i(i))n should be identified with the 
1— point velocity-frequency density function (VFDF). 
Namely, introducing the short- way notation f^f req ) (t) = 
fl^ req \ri, v,i, a), it should result 



(17) 



where 



f[ freq \t) = Jim 1 Y, Nx^^a-Z) (18) 

N— >oo iv L — ' 



i=l,N 



and Nx(ri, v,t, a; Z) is the frequency associated to the 
fluid velocity field V(r,i, a) [see related discussion in Ap- 
pendix B and in particular Eqs. (1114[) - (|118[) ]. Hence, by 
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definition, from Eq. (fT"7f it follows that it must be 

jd 3 v(f 1 (t)) n = J d^f\ freq) it) = l. (19) 
u u 



lC.b - The deterministic limit of /i 

Regarding the search of a deterministic limit for /i, 
there is manifestly only one case when this can happen 
and it occurs in the limiting case of the stationary solu- 
tion (of INSE) Z a = {V = 0,pi = 0} (or null solution), 
for which f\ = . In fact, in such a case f\ manifestly 

must coincide with a distribution, i.e., 

lim /i = *(v) (20) 

|V|->0+ 

(deterministic limit), while also there manifestly results 
lim fl freq \t) = <5(v). (21) 

|V|->0+ 

Since the null solution can always in principle be 
reached (i.e., moving backward in time and by applica- 
tion of a combination of appropriate volume forces and 
boundary conditions acting on the fluid, both to be suit- 
ably defined), it follows that Eq. (f20"|) might always be 
regarded as a possible alternative initial condition for f\ . 
In the following we intend to prove that, as an example, 
it can be represented as the limit function 

*(v) = lim /m(v;pi) (22) 

|V|->0 + 

of a Gaussian PDF 

/M(v-V;pi) = ^^expj- l|v ~ V|1 1, (23) 

^ v ih, P { v th P J 

1/2 

where V = and Vthp — (2pi/p ) is the thermal ve- 
locity due to pi [see Eq.fl27}]. It is obvious, however, that 
the previous requirement does not provide a unique func- 
tional form for f\. For instance, when both V(r, t, a) 
and p±(r,t,a) are considered infinitesimals (of order e), 
Eq. ([20|) only requires that 

h(t) = _Mu;pi,or) + 5fi(r,u,t,a; Z), (24) 

with Sfi infinitesimal of order 0(e). 

In the following {/i,r} will be denoted as statistically 
complete if the PDF fx : 

A) admits for all (r,t) <E fl x I (including the initial 
time t ) and G — 1, v,u 2 /2, uu, uu 2 /2 the velocity and 



phase-space moments J dvGfi and J dvfi In fi and sat- 

u r 
isfies the constraint equations 

JdvGh = l,V(r,t,a), Pl (r,t,a), (25) 

u 

S(fi(t)) - St. (26) 

Here pi(r,t, a) > denotes the kinetic pressure 

Pi{r,t,a) =p(r,t,a) +p (t,a) -<p(r,t,a), (27) 

withpo(^i ot) > (the pseudo-pressure) a strictly positive, 
smooth, real function and </>(r, t, a) a suitably defined 
potential; 

B) at the initial time t = t D satisfies the constraint 

(TO; 

C) satisfies the constraint defined by the deterministic 
limit flH|). 

The constraint equations ([2l ty -([?? ]) .(|r7] ) and ^ are 
here denoted as physical realizability conditions. 



ID - Open issues 

Here we shall consider a class of statistical models 
{/i,r} , based on the 1— point PDF, fi [with fi defined 
on r and T identified with ([3"])]. which yield a complete 
inverse kinetic theory for the INSE problem. In other 
words, each {/i,r} should yield the complete set of fluid 
fields, {Z}, to be represented in terms of suitable veloc- 
ity moments of the same PDF. In addition, {/i,r} will 
be required to hold for arbitrary fluid fields {Z} which, 
in the domain of existence Q x I, are strong solutions 
of the INSE problem, the latter - for greater generality 
- to be considered either deterministic or stochastic (see 
Appendix A). 

In the construction of statistical models of this type 
several interesting issues arise, which are related both to 
the prescription of the initial conditions on the PDF and 
to its time evolution. In particular, they concern whether 
there exists a statistical model {/i, T} , such that: 

1. (Problem 1: uniqueness condition) both {/i,r} 
and the PDF fi are unique; 

2. (Problem 2: conditions of statistical completeness) 
{/i,r} is statistically complete, namely f\ satisfies 
the requirements posed in Subsection 1C. 

Here we shall consider, in particular, the case in 
which the frequency fi^ req \t ) is an ordinary func- 
tion (i.e., not a distribution), consistent with the re- 
quirement that /i is an ordinary function too. The 
second requirement is that in the limit in which 
pi -> 0+ and |V| -> 0+, there results (|2U)) . 
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3. (Problem 3: moment- closure condition) can be de- 
fined in such a way that it satisfies a moment- 
closure condition. In other words: whether there 
exists a finite set of moment equations of IKE © 
which are closed. This is actually a basic require- 
ment of IKT. Hence, it should be satisfied by con- 
struction, if {/i,T} relies on IKT pj . 

4. (Problem 4-' kinetic-closure condition) it can be de- 
fined in such a way that at any time t > t (with 
t <E I), the statistical equation advancing in time f\ 
depends, besides f\ , only on a finite number of mo- 
ments of the same PDF , which include necessarily 
the complete set of fluid fields {Z(r,t)} . In such a 
case {/i,r} is said to satisfy a kinetic closure con- 
dition. This assumption is - in some sense - analo- 
gous to the closure problem for the ML hierarchy. 
In both cases it effectively involves the construction 
of the mean-field force which advances in time the 
1-point PDF. 

5. (Problem 5: determination of multi-point PDF's) 
the statistical model {/i,r} can be constructed in 
such a way that determines uniquely the multi-point 
PDF's, as well as the related observables. 

6. (Problem 6: closure condition of the statistical 
equations for multi-point PDF's) the statistical 
equations for the multi-point PDF's depend only 
on fx. 

In the remainder a statistical model which satisfies 
both the moment and kinetic closure conditions indicated 
above in Problems 3 and 4 will be denoted as closed. 



IE - Goals and scheme of presentation 

Here we claim that the IKT statistical model can be 
defined in such a way to provide an explicit solution of 
problems 1-6. The construction of the 1— point PDF 
is achieved adopting the so-called inverse kinetic theory 
(IKT) for fluid dynamics (Tessarotto et aL_2004-2009 
0> [l^ [l^ [l^ l2M> [2^ [IB, IH^, [3^, |3^| ). 
This type of approach can be formulated both when the 
fluid is considered regular and turbulent, i.e. the corre- 
sponding fluid equations (INSE) and the fluid fields {Z} 
are respectively considered deterministic and stochastic 
[lH, In particular, f\ is defined in such a way that: 

A) both the PDF and the related mean-field force (F) 
which determines its time evolution [see Eq.(|8|)] are con- 
ditional observables. This concept [see Appendix A, Sub- 
section A.l] is shown to imply the uniqueness and closure 
properties of the statistical model; 

B) the PDF satisfies suitable physical realizability con- 
ditions. This requires, in particular, that the complete 
set of fluid fields {Z} is necessarily represented in terms 



of moments the PDF [ll|, |lj, [lj , while the initial PDF 
must be defined so that its configuration-space average 
is suitably prescribed [see Eq. (Tl7|) ]. In such a case, the 
solution of the initial condition for fi (see Problem 2) is 
uniquely achieved by invoking the Principle of Entropy 
Maximization (PEM; Jaynes, 1957 [11]). This permits to 
identify two possible solutions for the initial PDF, corre- 
sponding respectively to the initial condition 1a and to 
1a and 1b jointly. 

The first one is realized by a local Gaussian PDF (/m)- 

The formulation of the corresponding statistical model 
{/Af,r} and the analysis of its basic properties are dis- 
cussed in Section 2. First, it is pointed out that PEM re- 
quires necessarily that f\ = fM must be a conditional ob- 
servable for all t <E L Then, by identifying the mean-field 
force F(/a/) with a conditional observable, it is proven 
that it is uniquely defined, both as a function of the 
kinetic velocity v (actually a polynomial of second de- 
gree with respect the relative kinetic velocity u = v V) 
and of the fluid fields {Z} = {V,pi} . As a consequence, 
the statistical model {/m, T} is unique (THM.l; see also 
Problem 1). Moreover, the complete set of fluid fields 
{Z} are uniquely determined as moments of f^j, which 
means that the classical dynamical system defined by the 
initial-value problem (|7|) determines uniquely both the 
time-evolution of the PDF and of the complete set of fluid 
fields {Z}. However, unless the initial constraint (|17| is 
fulfilled by fm (which is generally not the case), {/mjT} 
is not statistically complete (see Corollary to THM.l). 

The construction of the statistical model P} which 
takes into account such a case is carried out in Section. 
3. The condition of statistical completeness requires 
that the initial PDF satisfied simultaneously conditions 
1a and 1b- It is found that this generally provides (at 
t = t a ) a non-Gaussian initial PDF of prescribed form 
(see THM.2). Its time evolution depends now on the 
mean-field force F(/i). Under the assumptions that its 
dependence in terms of the kinetic velocity v can only 
occur, as in the Gaussian PDF, only via a polynomial 
of second degree in u and identifying F(/i), as in the 
previous case, with a conditional observable its form is 
found to be unique as the statistical model T} while 
{/i,r} is also statistically complete. 

In the same section (Subsection 3D) the closure prob- 
lem (Problems 3 and 4), based on the IKT statistical 
model, is formulated. This refers, in particular, to the 
construction of its formal solution (i.e., see also para- 
graphs IB and ID). This is achieved, both for determin- 
istic and stochastic fluids described by the INSE problem 
(as defined in Appendix A). In our theory this is done 
by determining directly f\ , rather than its stochastic 
average (/i) . As a basic consequence, it follows that the 
statistical equation for fx is necessarily closed, namely it 
depends only on f\ and a finite number of moments of 
the same PDF. The proof of the closure property of the 
IKT statistical model is proven in THM.3. 
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There are several new contributions and basic conse- 
quences of the theory here presented. 

Besides the (generally non-Gaussian) characterization 
of the initial PDF, it is found that, both for regular and 
turbulent flows, its time evolution is provided by a Liou- 
ville equation, while the corresponding Liouville operator 
depends only on a finite number of velocity moments of 
the same PDF. Hence, its time evolution depends (func- 
tionally) solely on the same (1-point) PDF. An interest- 
ing issue is also provided by the comparison, carried out 
in section 4, between {/i,r} and the common statisti- 
cal model, denoted {/ij,r}, laying at the core of the 
customarystatistical approaches [i.e., the HRE [J d, @| 
and ML [7J, |8| approaches]. The latter, although both 
unique and closed (in the sense of Problems 3 and 4), is 
proven to be statistically incomplete (see THM.4), since 
the corresponding PDF cannot generally fulfill the phys- 
ical constraint placed on it by Eq. (fT7|) [at t = t Q or at 
any time t 6 /]. 

The connection with previous statistical approaches 
is investigated. In particular, it is shown that f\ can 
be identified with a suitable stochastic average of the 
PDF fu (see THM.5). A remarkable consequence of 
the present theory is that it affords the exact construc- 
tion of the corresponding statistical equation for the 
stochastic-averaged PDF. (Section 4). The latter is shown 
to depart from the customary transport equation consid- 
ered in the literature [see, for example, Dopazo [Bj] and 
Pope [g]]. The relationship with the Hopf functional- 
differential method [l[ is also displayed. 

Finally, the explicit construction of multi-point PDF's 
and of the related observables (Section 5) is achieved. 
In particular, the statistical evolution equations for the 
multi-point PFD's are shown to maintain the form of 
Liouville equations. As a practical application, the ex- 
plicit construction of reduced 2— point PDF's - and of 
their related statistical equations - both usually inves- 
tigated in experimental/numerical research in fluid dy- 
namics, is presented. In fact, despite not being them- 
selves observables, they are nevertheless related to phys- 
ical observables (or conditional observables), such as the 
velocity difference between different fluid elements, usu- 
ally adopted for the statistical analysis of turbulent flu- 
ids. Finally, in Section 6 the conclusions are drawn. 



2 - IKT APPROACHES 

A fundamental aspect of fluid dynamics is the con- 
struction of statistical models {/i,r} in which the 
1— point PDF is the solution of a so-called inverse prob- 
lem, involving the search of a so-called inverse kinetic the- 
ory (IKT) able to yield the complete set of fluid equations 
for the fluid fields. A particular realization for {/i,r} 
is provided by (the already mentioned) ML approach, 



which is based on the position (jlip and the particular 
solution ([121 . In such a case it follows, by construction, 
that f\ depends explicitly, and not just merely in a func- 
tional sense, on the fluid field V(r,<). Hence, IKE © 
implies necessarily INSE (and therefore can be viewed as 
an inverse kinetic equation). Nevertheless, it is obvious 
that the fluid pressure p(r,t) cannot be represented as a 
moment of the same PDF. 

In this connection, however, a more general viewpoint 
is represented by the search of so-called complete IKT's 
able to yield as moments of the PDF the whole set of 
fluid fields {Z} which determine the fluid state and in 
which the same PDF satisfies a Liouville equation. This 
implies that in such a case there must exist a classical 
dynamical system, of the type defined by Eq. fiSp, whose 
evolution operator Tt ot t advances in time both the PDF 
and the related fluid fields, while preserving - at the same 
time - a suitable probability measure (Frisch, 1995 |19|). 
Despite previous attempts (Vishik and Fursikov, 1988 
[20l | and Ruelle, 1989 [5l[) the existence of such a dy- 
namical system has remained for a long time an unsolved 
problem. 

This type of approach has actually been achieved for 
incompressible NS fluids [HI, [13], with the discovery 
of the (corresponding) NS dynamical system which ad- 
vances in time the complete set of fluid fields {Z} . Its ap- 
plications and extensions are wide-ranging and concern 
in particular: incompressible thermofluids (25[, quan- 
tum hydrodynamic equations (see [23], [26|), phase-space 
Lagrangian dynamics 27j, tracer-particle dynamics for 



thermofluids 28|, 32|, the evolution of the fluid pressure 
in incompressible fluids [29j ] , turbulence theory in Navier- 
Stokes fluids [l3|, 25] and magnetofluids [3] and appli- 
cations of IKT to lattice-Boltzmann methods [3l|. In 
the following we intend to investigate, in particular, its 
consequences for the problems posed in this paper. 



2A - The IKT statistical model - Basic assumptions 

Let us now show how a statistical model of this type 
for the 1-point PDF [i.e., {/i,T}], which fulfills the re- 
quirements posed in Problems 1-5 and holds both for 
regular and stochastic flows , can be achieved by 

suitably modifying the IKT approach earlier developed 
by Tessarotto and coworkers [Hi [13] (see also Ref.pjl). 

Such a theory, it must be stressed, is based on some 
of the axioms (such as the conservation of entropy, the 
principle of entropy maximization or the regularity as- 
sumptions) which are typical of CSM. In particular, the 
IKT of Refs. [13, E3]) requires that fx(t) = f x {y.,t,a;Z) 
is a particular solution of the 1— point IKE J9]) and satis- 
fies the following assumptions (Axioms #0-#4), impos- 
ing that for all (x,<) eT x I : 



(Axiom #0: regularity) fi(t) is an ordinary, 
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strictly positive function. By assumption it is mea- 
surable, i.e., it admits the velocity- and phase- 
space moments 

JdvGfi (28) 
u 

respectively for G = 1, v, p Q u 2 /3, uu, mt 2 /3 and for 

8(fi(t)) = -/djc/i(r,u,t > a;Z)ln/i(r,u,t,a;Z), 
r 

the so-called Boltzmann-Shannon (BS) entropy. 
Moreover, fx is suitably smooth, in the sense that it 
is is continuous T x /, is a particular solution of IKE 
© and its velocity moments G 2 ,Gz = v,p u 2 /3 
are respectively of class 



• (Assumption #1 : conditional observables) f\ and 
the mean-field force F are conditional observables; 

• (Assumption #2: statistical completeness) {fi,T} 
is statistically complete, namely fi, satisfies the 
physical realizability conditions placed by: (As- 
sumption #2a) the velocity moments (|2l)j) -(f2"T |) ; 
(Assumption #2b) the VFDF, i.e., Eq.|Q2|>; (As- 
sumption #2c) the deterministic limit (|2T))) . 

In the remainder Axioms #0-4 will be assumed to hold, 
together with Assumption #1 and #2, for the statistical 
model {/i, r}. 

2B - The Gaussian particular solution 



G 2 eC^ 1 Hn x I), 
G 3 eC( 2 '°)(f! x /); 



(29) 



• (Axiom #1: principle of correspondence) fi(t) sat- 
isfies identically in x I the constraints (j2"5j) - (j2"T|) 
and the fluid fields (including the kinetic pressure 
pi) can be considered (conditional) observables. As 
a part of the same axiom it is required, further- 
more, that /i(r, u,i; Z) satisfies suitable kinetic ini- 
tial and boundary conditions [see Ref.]. The latter 
are defined so that, by construction, the moments 
(f2"5")) - (j2"7T) identically satisfy the initial and bound- 
ary conditions prescribed by the INSE problem [see 
Eas.(|9"2j) and ([53]) in Appendix A]; 

• (Axiom #2: moment- closure condition) the mo- 
ments equations, corresponding to G — 1, v, p Q u 2 /3 
and evaluated in terms of IKE ([5]), define a system 
of closed equations, which coincide with the com- 
plete set of fluid equations provided by INSE [see 
Eqs.{88])-{90j in Appendix A]; 

• (Axiom #3: principle of conservation of entropy, 
or constant H-theorem) f\(t) satisfies the con- 
straint equation [also known as constant H-theorem 



(30) 



• (Axiom #4-' maximum entropy principle) f\(t) sat- 
isfies, at t = t D , the constrained maximal variational 
principle (also known as principle of entropy max- 



imization or PEM; Jaynes, 1957 
5S(h(t)) = 0. 



33) 



(31) 



In this paper we shall impose, furthermore, two new 
assumptions introduced to satisfy the previous problems 
(Problem 1-7), requiring that: 



Let us now show that imposing PEM (at t — t a ) and 
requiring solely the specification of its functional class 
{/i} - i.e., not imposing also the validity of the initial 
constraint equation (fl7|) for the 1-point PDF - uniquely 
determines its initial value fi(t ) = f(y.,t ,a;Z). 

For definitcness, let us assume that fi(t ) is an ordi- 
nary, strictly positive function, requiring initially that it 
satisfies only the constraint provided by the fluid fields, 
namely Eq. (|9"2"|) [Assumption #2a] and by the determin- 
istic limit (|2"0|) [Assumption #2c]. Invoking the method 
of Lagrange multipliers, the PEM variational principle 
[ (|3ip ] implies that at the initial time t a and for arbitrary 
variations 5fi(r, u,i, a; Z) it must result identically: 

J dxSf 1 (r,u,t,a;Z){l + lnf 1 (r,u,t,a;Z)+ (32) 
r 

+A G + Ai -u + A 2 u 2 } = 0. 

Here A Q ,Ai and X 2 are Lagrange multipliers to be de- 
termined by imposing the correspondence principle [i.e., 
the moment equations (|25p]. It follows that: 1) at 
t — t , /i(r, u,t, a; Z) necessarily coincides with a Gaus- 
sian distribution (|23|) [l2l [23| carrying the fluid veloc- 
ity V(r,i D , a) and the kinetic pressure pi(r,t a , a); 2) Ax- 
iom #0 implies that pi(r, t ) must be strictly positive; 
3) /m(^o) is necessarily a conditional observable, so that 
when a and v (or the relative kinetic velocity u) are con- 
sidered prescribed, /m (u = v V(r,i D , a);pi(r, t D ,a)) is 
unique and depends in a prescribed way only on observ- 
ables; 4) moreover, due to the arbitrariness of the choice 
of the initial time t a , the Gaussian PDF (|2"3"|) is necessar- 
ily a particular solution of IKE © for all t 6 / ; 5) as 
a consequence, also the same equation at all times is a 
conditional observable in the sense indicated above. 

On the other hand, imposing for Eq.Q the particular 
solution 

/i(r,u,t,a;Z) = f M (u;px(r, t, a)) (33) 

implies that the vector field F must depend functionally 
on f M and {Z} , i.e., it is of the type F(f M ) = F(x,£; f M ) 
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|12l . |23| |. One finds, however, that the most general 
admissible form of the vector field F(/m) [namely one 
which is consistent with the requirement ([55]) ]. is of the 
type: 

F(x,t, a; f M ) = F (x,t, a; f M )+ (34) 
+Fi(x,t,a;/ M ) + AF(x,i,a;/ M ). 

Here AF is an arbitrary "gauge" vector field satisfying 
the homogeneous equation 



^ • (AF/m) = 0, 



(35) 



i.e., it does not contribute to IKE when Eq. (|33|) holds 
identically, while F , Fj and A(r,i; Jm) are respectively 
the vector and scalar fields 



F (x,t,a;/ M ) = — fe+ 
+i u • vv + ivv • u+^V 2 V, 



Fi(x,i,a;/ M ) = -uA(r,t;f M ) 



(36) 



(37) 



-th 



Vlnpi 



-4(r,t,Q!;/. 



: 7/i 



i a 



A/ J 



Pi Ot pi 
1 



y 2 / 2 + (38) 



+V • VV A /2 V • f - z/V-V 2 V 

Po 

Then, by construction Fo and Fi are conditional observ- 
ables. In particular Eq. ([35|l requires 



(39) 



AF(x,i, a; f M ) = u • E = AFi(x,i, a; f M ), 



where E = E(r,t) is an arbitrary antisymmetric second- 
order tensor. From Eg. (|59"|) it follows that AF is mani- 
festly non-observable. In fact, introducing the transfor- 
mation 



AF ^AF'=fcAF 



(40) 



with k G K arbitrary and non-vanishing, it yields for F 
an admissible form too [i.e., consistent with Eq. ([33]) ]. 



2C - Properties of {/m,T} 

Let us now pose the problem of resolving the indeter- 
minacy of F(/m). In Refs. [H, 0| for Gaussian solu- 
tions the uniqueness of F was achieved based on the re- 
quirement of consistency with extended thermodynam- 
ics, namely by imposing a suitably-prescribed form for 



higher-order moments of the Liouville equation. In par- 
ticular this was found to require 



AF = 0. 



(41) 



Here we point out, however, that Eq. (|41[) is equiva- 
lent to impose that F(x,i, a; /m) is a conditional observ- 
able (see Assumption #2), namely a uniquely-prescribed 
(polynomial) function of the kinetic velocity v (as well of 
the variables and r,t, a). Hence, the following theorem 
holds: 

THM.l - Uniqueness of {f M ,F} 

In validity of Assumptions #1, 2 and 3a, 3c the statis- 
tical model {/mjT} defined by Eas. &23\) , with mean-field 
force F(x,i, q;;/m) prescribed by Eqs.l[34\)- $38\) and sub- 
ject to the constraint |^_?[ ) ; is unique. 

PROOF 

The proof is immediate. Uniqueness follows, in fact, 
besides the axioms of IKT (Axioms #0-4), from the 
requirement that both /m and F(x,i, q;;/m) be condi- 
tional observables and hence, in particular, from Eq. (j4Tj) . 
Q.E.D. 

In Ref. [ijj the statistical model {/Af,r} was proven 
to determine uniquely, thanks to Axiom #1, the complete 
set of fluid fields. Such a result holds, however, in princi- 
ple for an arbitrary choice of AF of the form (|59")) , namely 



also in the case in which AF ^ 0. In Refs. [2J, [24 1 
the uniqueness of the mean-field force for the Gaussian 
PDF (/m) was achieved based on phenomenological argu- 
ments, i.e., the comparison with extended thermodynam- 
ics. The present result shows, however, that uniqueness 
tor F(x,t, a; /m), and hence {fm,F} too, can actually 
be achieved based on the physical prescription that both 
quantities are conditional observables (see Assumptions 
#1 and 2). 

Nevertheless, it is obvious that the statistical model 
{/m,T} generally is not statistically complete, in the 
sense indicated above. In fact (generally) /m does not 
satisfy the constraint imposed by the initial condition 
(fiTl There it follows immediately: 

COROLLARY to THM.l - Statistical incom- 
pleteness of {/m,! 1 } 

In validity of THM.l the statistical model {/Af,r} is 
generally statistically incomplete. 

PROOF 

We notice that for a Gaussian PDF (|23| the determin- 
istic limit ([20| exists. This requires letting p\ — > + and 
V| — » + , implying also that p\ must be uniquely deter- 
mined (i.e., it is necessarily an observable). However, it 
is obvious the constraint (|9"2"|) may not be fulfilled by /m- 
In fact, there results generally aX t — t a (for finite p\ and 

ivn 



fl fre "\v,t, a ;Z) 
(/A/)n 



^1- 



(42) 
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Q.E.D. 

This leaves open the question of extending IKT in such 
a way to satisfy this requirement (of statistical complete- 
ness) . In the remainder we intend to show that this con- 
straint generally implies a non-Gaussian initial PDF. 



3 - IKT STATISTICAL MODEL: NON-GAUSSIAN 
CASE 

An interesting question is posed by the case in which 
the PDF at t = t a , fi(t) satisfies also the initial condi- 
tion determined by Eq. |J7| ) [Assumption #2b] . If the 
fluid velocity V(r,i) is bounded in the domain fi, the 
constraint (|17| implies necessarily that the subdomain of 
velocity space in which /i remains strictly positive is gen- 
erally a bounded subset of M 3 . Hence, this prescription 
generally corresponds to an initial 1-point PDF which is 
locally non- Gaussian, specifically because of the (miss- 
ing) tails of the PDF. 



3A - Solution of the initial-value problem for /i 

In this case, let us consider (at t = t Q ) fi of the general 
form 



2. unless there results identically 



(/On 



1, the initial 



h{t) = </i(*)> f 



h(t) 



(43) 



Ht))n 

with (/i(i ))n determined by Eq. (fl~7"l) and h(t D ) to be as- 
sumed again as a strictly positive and regular real func- 
tion. Then considering variations of /i the type 

i.e., defined so that there results identically <5(/i) n = 
S (h) n = 0, the variational principle (|3~lj) (PEM) requires 
in this case that at the initial time t Q , h(t a ) must fulfill 
the variational equation 

f dx5h(t )^^{l + \nh(t )+ (45) 
r 

+A + Aj -u + A 2 u 2 } =0. 

Thus, /i (t Q ) is necessarily of the form (|4"3")) , while at t = t Q 
the function h(t) reads 

h(t) = cxp{-l - A - Ai • u-A 2 u 2 } , (46) 

with the Lagrange multipliers A D ,Ai and A 2 to be de- 
termined again imposing the moment equations (|25p . It 
follows that: 

1. due to the arbitrariness of the choice of t a , thanks 
to Axioms #3 (conservation of the BS entropy) and 
4 (PEM), it follows that for all t e /, fi(t) is nec- 
essarily of the form (I43[) : 



PDF fi(t ) is generally not a Gaussian, and hence 
of the form fi(t) = gi(£)/M(i),with gi(t) ^ 1 to be 
assumed a suitably smooth ordinary function; 

3- fi(t ) is necessarily a conditional observable (As- 
sumption #1). 



3B - The Liouville evolution equation for /i 

To determine the time evolution of the PDF we require 
again that the mean field force F(/i) = F(x,t, a; fi) is 
a conditional observable. To determine explicitly F(/i), 
lest us require, first, that when Eq.([33|) holds identically, 
F(/i) must coincide with F(/ M ) [defined by Eqs. (|34"|) - 
(|39|) ]. The same manifestly must occur if f\ remains, in 
the whole phase space T, suitably "near" to Jm- In fact, 
in this case in the same set F(/i) and F(/m) must re- 
main close too (again in some suitable asymptotic sense). 
This delivers for F(/i) a prescribed polynomial repre- 
sentation in terms of the relative kinetic velocity u. In 
principle, in fact, F(/i) might include also polynomials 
of higher-degree in u. These terms, however, must nec- 
essarily vanish identically for f\ = Jm and cannot con- 
tribute to the moment equations [of IKE] which yield 
INSE. Since their form remains therefore arbitrary, they 
are manifestly non-observables. Hence such terms are 
ruled out by the requirement of F(/i) being a conditional 
observable. This implies that F(/i) must be necessarily a 
polynomial of second degree in the relative kinetic veloc- 
ity u, whose precise form is prescribed by imposing that 
the moment equations of IKE for G(v, r, t) — 1, v, p Q u 2 /3 
must necessarily coincide with INSE [ Eqs. (f8"8")) - (f9"0")) in 
Appendix A] . Let us introduce now the mean-field force: 

F(x,*,a;/i) = F (x,i,a;/i) + F 1 (x,i,a;/i) + (47) 
+AF(x,t,a;/ 1 ), (48) 

AF(x,i, a;/i) denoting here an arbitrary gauge vector 
field of the form 

AF(x,t, a; f x ) = AF 1 (x,t, a; f M ) + AP 2 (x,t, a; A), 

defined so that: a) if f\ = fM, it coincides with 

AF( x> i,a;/i)=AFi(x,*,a;/ M ); 

b) it does not contribute to the moment equations of 
IKE, evaluated for G(v,r,t) = l,v,p Q u 2 /3, namely it is 
defined so that there results identically 



(49) 



Hence AF(x,t, a; fi) satisfies manifestly the gauge condi- 
tion (|4H|) . The indeterminacy of F(x,i, a; f\) can avoided 
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again by imposing again the requirement that it is a con- 
ditional observable (Assumption #2), namely 



AF(x,t,a;/i) = 0. 



Here the vector fields Fo and Fi are [12 1 



(50) 



(51) 



P (x,t, a; fx) = — [V-n - Vpi + f B ] 
Po — 

+i u • vv + ivv ■ u+^V 2 V, 



Fi(x,i, a; /i) = i u ( A(r,t; /)+-V • Q- (52) 
2 I Pi 



1 



vn -q + v-n - - 



w 3 



Pi ^ — j 2pi 2 
with Q and II denoting the moments 



Q =p fd 3 vu^h, 
II = po/ d 3 wuu/i, 



and the scalar field A(r,t, a; /i) defined by Eq. 



(53) 



3C - Properties of {/i,T} 

Let us now prove that the IKT statistical model {fx, V} 
satisfies both Problems #1 and 2 also for non-Gaussian 
PDF's of the type (|4"3"1) . As a consequence, the following 
theorem holds: 

THM.2 - Uniqueness and statistical complete- 
ness of {/i,r} 

In validity of Assumptions #1, 2 and 3 (i.e., 3a-3c), let 
us impose that the statistical model {/i,r} is defined by 
requiring that fx is a strictly positive ordinary function 
which satisfies the initial condition |^3| ) with |^6j ) and 
IKE 0), with the corresponding dean field force defined 
by Eqs. {4^ -^52 ]) . Then {fx, T} is unique and statistically 
complete. 

PROOF 

In fact, first, the initial condition (|4"3")) determined im- 
posing PEM together with the constraints (|92p and (JTTJ) 
is manifestly unique and hence a fi(t ) is necessarily con- 
ditional observable. Furthermore, thanks to Assumption 
#2 F(x,t; /i) is necessarily of the form provided by Eqs. 
([4"T)) - (f5"2"j) . i.e., unique, so that {/i, T} is unique too. Sec- 
ond, the initial PDF satisfies by construction also to the 
requirement posed by Eq. (|17[) . Let us now prove that 
fx fulfills also the deterministic limit ([20]) . This follows 
invoking Eq. (|2ip and noting that in the limit pi — > + 
and |V| — > + there results manifestly 



h(t) 
P1 Zo+ (h(t)) n 

|V|-Q+ 



lim 



1. 



(54) 



This proves that is statistically complete too. 

Q.E.D. 

We remark that: 

• the proof that {/i,r} provides a complete IKT for 
the INSE problem was reached previously in Ref. 
(l2l |. This follows by noting that, by construction, 
(thanks to Axioms #0-2) the fluid fields {V,pi} 
necessarily satisfy the INSE initial-boundary value 
problem [see Appendix A] , when they are identified 
with the two velocity moments of the PDF evalu- 
ated for G = v, p a u 2 /3. 

• THM.2 yields a solution to Problem #1, namely 
that the initial the 1-point PDF coincides with 
the observable defined by the initial VDFD 
f} freq \v 3 ,t ,a;Z) [seeEq.CIJ)]; 

• In particular, THM.2 (as also THM.l) holds both 
in the case in which the fluid fields are deterministic 
and stochastic, i.e., both for the deterministic and 
stochastic INSE problems [defined in Appendix A] ; 

• In case of THM.2 the assumption that both the 
1— point PDF fx and F are conditional observables 
(Assumptions #1 and 2) is ultimately demanded 
also by consistency with the physical requirement 
posed by the constraint (TTT)) . In fact, it is obvi- 
ous that otherwise the uniqueness of {/i,T} might 
not be warranted. In particular, this means that 
configuration-space average of the PDF, i.e., either 
the continuous or discrete averages (/i)q and fx 
might not be unique, and hence these quantities 
would not be observables; 

• In Refs. [23|, [24| the uniqueness of the mean- 
field force F(x,i, a; fx) for a non-Gaussian fx, was 
achieved again based both on phenomenological 
arguments, i.e., besides the comparison with ex- 
tended thermodynamics, the requirement that F 
depends only on the minimum number of higher- 
order moments of fx, defined so that: a) they pro- 
vide the correct fluid equations; b) they vanish 
identically in the case (f33|) . THM.2 shows that 
the uniqueness [of F(x,t, a; fx)] is achieved simply 
based on the physical prescription that both fx and 
F(x,t, a; fx) are conditional observables (see again 
Assumptions #1 and 2). 



3D - Solution of the closure problem 

As indicated above, a desired property of statistical 
models in fluid dynamics would be the (possible) fulfill- 
ment of suitable closure conditions, permitting to assure 
that the relevant PDF advances in time by means of a sta- 
tistical equation which depends solely on the same PDF. 
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Here we wish to investigate the closure problem earlier 
posed [Subsections IB and ID] and, in particular, that 
specified by Problems 4 and 5. It is immediate to prove 
that a formal exact solution for such closure problems is 
provided by the IKT statistical model {fi,T} defined in 
the previous sections 2 and 3. In other words, {/i,T} is 
necessarily closed, i.e., the following result holds: 

THM.3 - Closure properties of {/i,T} (Closure 
Theorem) 

In validity of Assumptions #1, 2 and 3 (i.e., 3a-3c), 
the statistical model {/i,r} , defined by imposing the ini- 
tial condition |^ j| ) with (JBG , and with /i(r, v,i) assumed 
as a strictly positive ordinary function satisfying IKE Q), 
with mean field force defined by Eqs. J^Tp - $52\) . is closed. 

PROOF 

Due to the prescription (|34|) - ([38| of the mean-field 
force F(x,i, a; f\) it is immediate to prove that the mo- 
ment equations of IKE for G = 1, v, p Q u 2 /3, namely 

J d 3 vGL/i = (55) 
u 

coincide respectively with Eq. (l89|) (for the moments G = 
1 and p u 2 /3) and with Eq. ([90|) (for the second mo- 
ment G = v), i.e., with the complete set of PDE de- 
fined by INSE [see Eqs. flHHJ}- [MD, in Appendix A]. The 
latter, by assumption, define a closed system of equa- 
tions, hence the moment-closure condition (Problem 4) 
is necessarily satisfied. Furthermore, to prove that also 
the kinetic-closure condition (posed by Problem 5) holds 
it is sufficient to notice that the same mean-field force 
F(x,i, a; f\) depends, by construction, only on the fluid 
fields V,pi, Q,and n. Therefore, also the requirement of 
kinetic-closure is necessarily satisfied. Q.E.D. 

4 - COMPARISONS WITH PREVIOUS 
APPROACHES 

Interesting issues are posed by the comparison with 
previous statistical treatments and in particular the sta- 
tistical model {fa, T} , underlying both the HRE 0,0,01 
and ML 0, [1| approaches. This is relevant in special ref- 
erence to: 

• analyze basic properties of {/ir,r} (see Subsection 
4A); 

• determine the explicit relationship between the 
PDF's fi and fjj, which characterize the two ap- 
proaches (Subsection 4B); 

• the construction of the statistical equation for the 
stochastic-average of f\ (Subsection 4C), 

• the comparison with the HRE functional- 
differential approach (Subsection 4D). 



4A - Properties of {fu,T} 

It is immediate to show that {fn, T} : 1) holds both for 
the deterministic and stochastic INSE problems [see Ap- 
pendix A]; 2) realizes an IKT for INSE, which is unique 
and closed; 3) it is (generally) not a complete IKT (in 
fact, manifestly, neither the fluid pressure, nor the ki- 
netic pressure, can be represented in terms of velocity 
moments of ///); 4) in addition, it is (generally) not sta- 
tistically complete. Indeed, since fn is a distribution, its 
configuration-space average cannot generally expected to 
agree with the the observable f\^ req \-v,t, a: Z) or VDFD 
(1— point velocity-frequency density function). The fol- 
lowing result holds: 

THM.4 - Statistical incompleteness of {/#,r} 
The statistical model {/#,r} does not generally fulfill 
the constraint J_?7| ). 
PROOF 

To reach the proof, let us evaluate the configuration- 
space average of fn- For definiteness, let us adopt the 
definition of discrete average provided by Eq. (fT6]) . Thus, 
by denoting Zi = {Vi(t),pi(t)} the average fluid fields 
in the i— th cell of O (and evaluated at position r^), it 
follows that the configuration-space average of fu reads 

(Mt)) a = jjr E <K V - V *(*))» ( 56 ) 

* i=l,N, 

i.e., (/ff)n (v, i) * s stiM a distribution. Hence, it cannot 
generally be identified with the observable f± (which 
by assumption here it is considered as an ordinary func- 
tion)! For comparison, instead, the complete IKT ap- 
proach yields instead: 

</i(*)>n = TF- E /*(*)• ( 57 ) 

* i=l,N. 

where in this case by definition f\ is an ordinary func- 
tion (and hence its configuration-space average (fi{t)) n 
can be identified with the observable/conditional observ- 
able /f (i), which - on the contrary is generally an 
ordinary function. Q.E.D. 

A similar proof can be achieved also adopting the con- 
tinuous operator (fT4")) . We remark here that: 

1. the setting based on the definition p^|) is actu- 
ally consistent with the physical measurement pro- 
cess of the corresponding 1-point velocity frequency 
[fi^ reQ \t)]-> based of the discretization of the fluid 
domain Q]; 

2. the extension of THM.4 to a generic stochastic 
model of the type {(//f),r} is not possible, as 
proven by the subsequent discussion [see in par- 
ticular THM.5]. 
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4B - Relationship between {/i,T} and {/ir,r} 

Let us now pose the problem of determining the re- 
lationship between the 1-point PDF which characterizes 
the IKT statistical model {/i,T} and the PDF f H (t) 
associated to the statistical model {fu,T}. We intend 
to show that fx [with fx = /i(r, u,i), u = v — V(r,i) 
and x = (r, v) G T] can be identified with a suitable 
stochastic-average of fjj, to be defined in terms of an ap- 
propriate nonhomogeneous and non-stationary stochastic 
model (see related definitions in Appendix A). 

To assess properly precisely the statement let us intro- 
duce for the fluid fields {Z} the stochastic representation 

{Z(AV)} = {V(r,i)+AV,p(r,i)}, (58) 

where {Z} = {V(r, t),p(r, t)} and AV G R 3 denote re- 
spectively an arbitrary particular solution of the INSE 
problem [see Appendix A] and an arbitrary velocity 
fluctuation. Here by assumption, the vector AV = 
(AVi, AV2, AV3) will be considered as a set of stochastic 
hidden variables characterized by the stochastic PDF 

5l (AV,r,i) = A(r,AV,£), (59) 

with fi denoting the 1-point PDF characterizing the 
IKT statistical model {/i,r} previously introduced 
[and obeying Eq. ([9]) when replacing A V — > v] , The set 
{AV,gi} defines therefore (a generally non-homogeneous 
and non-stationary) stochastic model defined so that, by 
assumption, its moments are necessarily defined so that 

(l) a = J d 3 AVA(r,AV,i), (60) 
= J d 3 AVAV/i(r, AV,i), (61) 

R 3 

\ Po v 2 th {v,t) = ^ |d 3 Avi(AV) 2 /i(r,AV,462) 

R 3 

Here the moment on the r.h.s. of Eq. (f6"2")) represents, 
up to the constant factor |p Q , the stochastic mean value 
of the stochastic kinetic energy per unit mass h (AV) 2 , 
while hpoVthfti ^) an< ^ v th(^, t) are the thermal energy and 
the corresponding thermal velocity produced by the kinetic 
pressure px (r, t) [with px (r, t) defined according to Eq. 

H27M 

/ 2pi(r, t) 

v th (r,t)=\ i . (63 

V Po 

It follows that: 

THM.5 - Representation of fx in terms of fu 

The stochastic representation i58\) and position A59\) 
can always be introduced. It follows that the stochastic 



average of fu, defined with respect to the stochastic av- 
eraging operator il06\) , reads identically 

(f H ) AV = fi(r,u,t), (64) 

with fu given by Eq.fiSty. 
PROOF 

First we notice that if {Z} is an arbitrary particular 
solution of INSE, also {Z(AV)} is manifestly a particular 
solution of the same equation. The proof is reached by 
introducing in the NS equation (j9"0)) [see Appendix A] a 
suitable stochastic volume force f + Af , with Af defined 
so that 

— Af = AV • VV. (65) 

Po 

Hence, the stochastic representation (|5"8"|) always holds, 
with AV to be considered as stochastic hidden vari- 
ables. In addition, since the definition of gi(AV,r,t) 
is arbitrary, the position (|59p can always be introduced. 
These definitions uniquely prescribe the stochastic model 
{AV,5i} and the related stochastic averaging opera- 
tor (| 1 06[) [see Appendix A]. It follows that (/h)av = 
J d 3 AV./i(r, AV,t)«5(v - V(r,t) - AV) yields Eq.©, 

R 3 

which is therefore identically satisfied. Q.E.D. 
In conclusion: 

• the 1-point PDF of the IKT statistical model can 
be considered simply as a possible stochastic real- 
ization of the PDF fu, achieved by means of the 
stochastic model {AV,gi} , in which g\ is properly 
related to the 1-point PDF characterizing the IKT 
statistical model; 

• in view of the positions (|62|) and (|63|) . {AV,<7i} can 

be viewed as the stochastic model which takes into 
account the thermal motion produced in a NS fluid 
by its kinetic pressure pi(r,t) [see Eq. (|27p ], 

4C - Consequences - Statistical equation for 

The same conclusion (i.e., THM.5) holds manifestly 
also in the case in which the fluid fields {Z}, as well 
as the same 1-point PDF fx and the mean-field force 
F(/i), are stochastic in the sense of Eq. (fl"3")) . i.e., they 
depend on a suitable set of stochastic hidden variables 
ol G V a C R™, here considered independent of AV and 
characterized by a stochastic PDF g. For definiteness, 
let us assume that g is homogeneous and stationary, i.e., 
of the form g =g{a). Then the ensemble average of the 
PDF, (/i) , can be identified with 

(fx) = (fi) a , (66) 

with fx defined by Eq. (f6"4"|) and the brackets (•) denoting 
the stochastic averaging operator (|106[) [Appendix A]. 
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Hence, as a consequence of THM.5 and IKE, it follows 
that (fi) obeys necessarily the stochastic-averaged IKE 

urn 



d_ 

dt 



(A) 



d 



<9v 



{(SF(Ji)Sfi)} 



{(F(/x)) </!>} = (67) 
(68) 



Here 6f± and (5F(/i) denote the stochastic fluctuations 
[see Eas. l|llip and pi2|) in Appendix A], while F(/i) is 
defined by Eqs. (|47 p ,(|50] ) .([CT ]) -(|55 ]) and §E§. We remark 
that Eq. (j6"T|) takes into account the constraints placed on 
by the physical realizability conditions (see Section 
2), as well the axioms of IKT (Section 2). Hence: 

• it is appropriate for describing homogenous and 
stationary turbulence in NS fluids; 

• however, it can be generalized, in principle, to 
a generally non-homogeneous and non-stationary 
stochastic PDF g, as appropriate for describing 
non-homogeneous and non-stationary turbulence 



The equation departs from the statistical (or "trans- 
port" ) equation [for (fi)] usually considered in the liter- 
ature [see, for example, Dopazo [H and Pope Q]. This 
fact is not surprising. In fact, unlike the customary ap- 
proaches, here: 

• the operator of ensemble average " () " , defined by 
Eqs. (|64|) and ([66]) . does not commute with the 
relevant differential operators (Jj, V, V 2 ); 

• by assumption, the 1-point PDF satisfies the Li- 
ouville statistical equation defined by ©, together 
with the physical realizability conditions imposed 
on the 1-point PDF [see Eqs. (J25j) - (J2TJ) , dTTJ) and 
([20]) ]. These constraints are not required in the 
customary approach 0, Q ; 

• the definitions of the vector field F(/i) here 
adopted satisfy the requirements placed by the ax- 
iomatic formulation (see Section 2, Subsection 2A). 



4D - Relationship between {/i,r} and the HRE 
functional-differential approach 



which implies 

f 

-c/)[y{x),t} = - dxy(x) 







v— A(x,t)+ (70) 



<9v 



{F(x,t;/i)/i(x,t)} 



(71) 



Then it follows that <p[y(x),t] must obey the following 
single linear functional-differential equation, manifestly 
equivalent to IKE (Liouville equation) (TT0)) . 



d f 

—(p[y(x),t\ = dxy(x)Q 



Sy(x) ' 



(72) 



where Q is related to the Liouville operator and is defined 

as 



Q- = v-|- + ^--{F(x ) t;/ 1 )-}. 
or ov 



(73) 



Eg. (f?2")) is analogous to the Hopf cj>— equation and has the 
explicit exact solution 



0[y(x),t} = <p[T to>t y(x ),t ] , 



(74) 



where T to j is the evolution operator associated to the 
Navier-Stokes dynamical system. 



5 - IKT FOR MULTI-POINT PDF'S 

The construction of multi-point PDF's is a problem of 
"practical" interest in experimental/numerical research 
in fluid dynamics. In fact, despite not being themselves 
observables, they are nevertheless related to physical ob- 
servables (or conditional observables), such as the ve- 
locity difference between different fluid elements usually 
adopted for the statistical analysis of turbulent fluids. 

In the present theory, unlike the ML approach, the 
statistical equation advancing in time the 1— point PDF 
fi [i.e., Eq.Q] satisfies, by definition, a kinetic closure 
condition. Hence, the construction of the multi-point 
PDF's is trivial. Nonetheless, it is still useful to analyze 
elementary implications (of the present theory) dealing 
with: a) the specific representation of certain "reduced" 
multi-point PDF's, defined in terms of the 1— point PDF; 
b) their dynamics, namely the statistical equations which 
they fulfill; c) their relationship with the relevant observ- 
ables. 



Finally let us consider the comparison with the HRE 
functional-differential approach [3, 0, 0]. In analogy to 
Ref. [HI let us introduce the functional 

<t>[y{x),t] = J djsyfxJACx.t), (69) 
r 



5 A - Liouville equations for the multi-point PDF's 

Let us assume, for definiteness, that fi(x.i,t, a; Z) is 
the 1— point PDF which is particular solution of the 
(1— point) Liouville equation Then, denoting fi(i) = 
/i(xi,t, a; Z) (for i = 1, s) the same PDF evaluated at the 
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states Xj = (rj, Vj) (for i = 1, s), the s— point PDF is the where V(i) = V(r,-,f), v 2 hp (i) — v 2 hp (r i: t) and de- 
probability density notes 



/ s (l,2,..s)= J] fi(i), (75) 

i— 1,5 

defined in the product phase-space V s = [] T. The sta- 

i=l,s 

tistical equation advancing in time f s follows trivially 
from Eq.©. In fact, denoting by F(i) = F(xi, t; fi) the 
mean- field force at the state x^ (for i = l,s) and intro- 
ducing the s— point Liouville operator (with summation 
understood on repeated indexes) 

L - (l,.. )a ) S * + v,A + A. {Fi(i)} , ( 76 ) 

it follows that / s (l, 2, ..s) satisfies identically the s— point 
Liouville equation (or IKE) 

L s (l,..,s)f s (l,2,..s) = 0. (77) 

5B - Consequences: reduced 2-point PFD's and 
2-point observables 

In terms of the 2-point PDF, f 2(1,2), a number of 
reduced probability densities can be defined in suitable 
subspaces of T 2 . To introduce them explicitly let us first 
introduce the transformation to the center of mass coor- 
dinates of the two point-particles with states (r^Vj) (for 
i=l,2) 

{r^Vi,!-^} -+ {r,R,v,V} (78) 

[with r = ^2^,R=^±^,v= and V =2a±a]. 

Then, these are respectively the local (in configuration 
space) velocity -difference 2-point PDF: 

Sa(n,r a ,v,t,a) =/ a d 3 V/ 2 (l,2) = 
= /d 3 V/ 1 (r 1 ,v + V,i, a ))/ 1 (r 2 ,V-v,f, a ;Z) ^ 

(defined in the space f2 2 x {/) and the velocity-difference 
2-point PDF, i.e., the observable 

%(r,v,t,a) = ( 52 (r + R,r-R,v,i,a)) Ri0 (80) 

(defined in T), (} R denoting the configuration-space av- 
erage operator acting on the center of mass. Hence, in 
terms of the average operator (| 1 4j) , there follows 

£(*.vA«) = ^//R (81) 
,g 2 (r + R,r-R,v,i, a). 

In the case of a Gaussian PDF [(021)], Eq.® delivers in 
particular the Gaussian PDF 



g 2 (ri,r 2 , v,t, a) = exp < § * 

(82) 



2 ^V P (1)+^V P ( 2 ) , eQ ^ 
«th = " 4 ■ ( 83 ) 

In a similar way it is possible to identify additional 
2-point observables. Precisely these can be defined as: 

1. the velocity- difference 2-point PDF for parallel ve- 
locity increments. Introducing the representations 
v = nw and r = nr, n denoting a unit vector, 
f2\\(i r , v,t) can be simply defined as the solid-angle 
average 

h\\(r,v,t,a) = J dSl(n)f 2 (r = nr,v = nv,t,a); (84) 

2. the velocity-difference 2-point PDF for perpendicu- 
lar velocity increments. Introducing, instead, the 
representations v = nv and r = n x br, n and b 
denoting two independent unit vectors, f2_i_(r,v,t) 
can be defined as the double-solid-angle average 

f 2± (r,v,t,a) = J dfi(n) J dfi(b) (85) 
/ 2 (r = n x br, v = nv,t, a). 

An interesting property which emerges from these 
definitions is that in all cases indicated above [i.e., 
Eqs.([80 ]) .(|84 )) and ([85]) ] the definition of g 2 given 
above [Eq. l(75j) ] implies that non-Gaussian features, 
respectively in / 2 ,/ 2 || and / 2 j_, may arise even if 
the 1— point PDF is Gaussian, i.e., the requirement 
(|33j) holds identically. This occurs due to velocity 
and pressure fluctuations occurring between differ- 
ent spatial positions ri and r 2 .More generally, how- 
ever, we can infer that - due to the constraint Q17p 
here imposed on the 1-point PDF - it is obvious 
that, if the fluid velocity V(r,t) is bounded in the 
domain f2, the same 1-point PDF, and hence the 
2-point PDF's, cannot be Gaussian distributions. 

5C - Statistical evolution equation for jfe 

From the 2— point IKE (177|) (obtained in the case 
s = 2) it is immediate to obtain the corresponding evolu- 
tion equation for the reduced PDS's indicated above. For 
example, the velocity-difference 2-point PDF / 2 satisfies 
the equation 

J o 
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It follows that, in particular, in the Gaussian case (|33[) 
this equation reduces to the (generally non-Markovian) 
Fokker-Planck equation 



d ? d 
or ov 



1 



dt dr J * dv jq 
F {T) g 2 (r + R, r R, v, t, a) = 0, 



(87) 



where the vector field F 



(T) 



pP(n : 



r 2 : V, t, a;/ M ) is 



given in Appendix C [see Eqs. (|129|) and ll30|) ]. 

An interesting issue is provided by the comparison 
with the statistical formulation developed by Peinke and 
coworkers 34, 35, 36, 37, [38J. Their approach, based 



on the statistical analysis of experimental observations, 
indicates that in case of stationary and homogeneous tur- 
bulence both the 2-point PDF's for parallel and velocity 
increments obey stationary Fokker-Planck equations. In 



particular, according to experimental evidence [37l . 
reasonable agreement with a Markovian approximation 
for Eq. I[87p - at least in some limited subset of parameter 
space- is suggested. Our theory suggests, however, that 
a breakdown of the Markovian assumption should be ex- 
pected due to non-local contributions appearing in the 
PDF's and in the corresponding statistical equations. 



6 - DISCUSSION AND CONCLUDING REMARKS 

An axiomatic approach, based on the IKT statistical 
model {/i, r}, has been developed for the statistics of the 
1-point PDF fi which characterizes an incompressible NS 
fluid. The paper contains several new aspects and basic 
consequences of interest in fluid dynamics and turbulence 
theory. 

Indeed, the theory here developed applies both to reg- 
ular and turbulent flows, characterized respectively by 
deterministic and stochastic fluid fields. In fact, in both 
cases the time evolution of fx is a Liouville equation 
(IKE) [see Eq.©] which evolves in time also the com- 
plete set of fluid fields (represented in terms of moments 
of the same PDF). 

In this paper an explicit solution of the problems 1-6 
posed in Subsection ID has been proposed. 

In particular, we have proven that - extending the sta- 
tistical approach earlier developed [ll|, \\2 - the IKT sta- 
tistical model {/i,r} can be uniquely determined. The 
present theory is based on two new hypotheses, i.e., A) 
that both fx and F(/i) are conditional observables and 
B) that fx satisfies suitable physical realizability condi- 
tions (see Subsection 2A). The first requirement permits 
to determine the mean-field force F(/i), while the second 
one uniquely prescribes - by means of PEM (i.e., Axiom 
#4) - the initial PDF fx(t ). In detail, we have shown 
that {fx, F} can be constructed in such a way to be: 

• unique (see THM.l, related to Problems 1); 



• statistically complete (see THM.2, Problem 2); 

• closed, i.e., it both moment-closure and kinetic- 
closure conditions (see THM.3, Problems 3 and 4) 

and furthermore: 

• that it determines uniquely all multi-point PDF's, 
as well as the related observables (see Problem 5); 

• that the statistical equations for multi-point PDF's 
depend only on fx and therefore, by definition, sat- 
isfy a closure condition (Problem 6). 

The theory has important consequences. 

First, it implies that the initial PDF is generally non- 
Gaussian (see again THM. 2). This conclusion holds even 
in the case in which the fluid fields are deterministic, 
namely for regular flows. In fact, the Gaussian 1-point 
PDF although unique (THM.l) does not generally pro- 
vide a statistical complete model {/m,F} (see corollary 
of THM.l). In addition: 

• thanks to the fluid and kinetic closure conditions 
imposed on the statistical equation for the 1-point 
PDF , i.e., IKE [Eq.©], fx depends only on a finite 
set of moments of the same PDF and its time evo- 
lution is independent of higher-order (multi-point) 
PDF's; 

• as a basic consequence, the exact statistical equa- 
tion for the ensemble-averaged (or stochastic- 
averaged) PDF (fx) has been obtained. This is 
found to be intrinsically different from the anal- 
ogous transport equation obtained in the past in 
the case of stationary and homogeneous turbulence 

The connection of the present theory both with previ- 
ous IKT approaches [ill Il2| and the HRE (Hopf, Rosen 
and Edwards [j], 0, [J) and ML (Monin and Lundgren 
fulfil) statistical treatments, has been pointed out (see, in 
particular, Section 4). Regarding, in particular, the last 
two approaches the following results have been reached: 

• the common statistical model, {/ij,r}, used in 
both approaches (HRE and ML) has been shown 
to be generally statistically incomplete (THM. 4) ; 

• the relationship between the fx and the PDF fn 
which characterizes the HRE and ML approaches 
has been determined. In particular, we have proven 
that fx can be identified with a suitable stochas- 
tic average of fn , via a generally non- homogeneous 
and non-stationary stochastic PDF [see Eq. ([64f ]; 

• the unique connection [via Eqs.(j64|) and (l66|) ] exist- 
ing between the ensemble-averaged PDF's (fx) and 
(fn) has been displayed; 
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• the relationship with the Hopf's functional- 
differential approach has been pointed out. 

Finally, as an application, explicit representations have 
been given for the reduced 2-point PDF's usually adopted 
for the statistical description of turbulent flows, rep- 
resented respectively by the velocity-difference 2-point 
PDF /2 and the velocity-difference 2-point PDF for par- 
allel and perpendicular velocity increments f2\\ and 
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APPENDIX A: THE MATHEMATICAL 
DESCRIPTION OF INCOMPRESSIBLE NS 
FLUIDS 



In fluid dynamics the state of an arbitrary fluid sys- 
tem is assumed to be defined everywhere in a suitable 
extended configuration domain ft x I [ft denoting the 
configuration space and ICR the time axis] by an ap- 
propriate set of suitably smooth functions {Z} , denoted 
as fluid fields, and by a well-posed set of PDE's, denoted 
as fluid equations, of which the former are solutions. The 
fluid fields are by assumption functions of the observables 
(r,t), with r and t spanning respectively the sets ft and /, 
namely smooth real functions. Therefore, they are also 
strong solutions of the fluid equations. In particular, this 
means that they are are required to be at least continu- 
ous in all points of the closed set ftx I, with ft = ft U dft 
closure of ft. In the remainder we shall require, for defi- 
niteness, that: 



1. ft (configuration domain) is a bounded subset of 
the Euclidean space E 3 on M 3 ; 

2. I (time axis) is identified, when appropriate, cither 
with a bounded interval, i.e., I=]to,ti[ C R, or 
with the real axis M; 

3. in the open set ftx I the functions {Z} , are as- 
sumed to be solutions of a closed set of fluid equa- 
tions. In the case of an incompressible Navier- 
Stokes fluid the fluid fields are {Z} = {V,p, S T } 



and their fluid equations 



P 

V- V 

iVV 

di ST 

Z(r,t ) 
Z(r,t)\ m 



Po, 

0, 
0. 

0. 

Zo(r), 
Z w (r,t)\gn ■ 



(89) 
(90) 

(91) 

(92) 
(93) 



where Eqs. (|88jl - (|9ip denote the incompressible 
Navier-Stokes equations (INSE) and Eqs. (|88p ~ 
(|93| the corresponding initial-boundary value INSE 
problem. In particular, Eqs. (f9"3"|) are respec- 

tively the incompressibility, isochoricity, Navier- 
Stokes and constant thermodynamic entropy equa- 
tions and the initial and Dirichlet boundary con- 
ditions for {Z} , with {Z (r)} and { Z w (r,t)\ gn } 
suitably prescribed initial and boundary-value fluid 
fields, defined respectively at the initial time t = t a 
and on the boundary dft. 

4. by assumption, these equations together with ap- 
propriate initial and boundary conditions are re- 
quired to define a well-posed problem with unique 
strong solution defined everywhere in ftx I. 



Here the notation as follows. 
operator 

D 



N is the NS nonlinear 



NV = — V-F H , 
Dt 



(94) 



with -gjV and Fh denoting respectively the Lagrangian 



fluid acceleration and the total force per unit mass 



D 
Dt 



V 



d_ 

dt 
1 

Po 



V- 



-Vp- 



V- vv, 
1 



Po 



-f + uV^V, 



(95) 
(96) 



while p > and v > are the constant mass density and 
the constant kinematic viscosity. In particular, f is the 
volume force density acting on the fluid, namely which is 
assumed of the form 



-V0(r, t) 



(97) 



</>(r, t) being a suitable scalar potential, so that the first 
two force terms [in Eg. ([TO]) ] can be represented as — Vp+f 
= — Vp r + in, with 



p r (v,t) =p(r,t) - 4>(r,t), 



(98) 



denoting the reduced fluid pressure. As a consequence of 
Eqs.|88]),(!89]) and (ffl) it follows that the fluid pressure 
necessarily satisfies the Poisson equation 

V 2 p = S, (99) 

where the source term S reads 

S = - Po V ■ (V • VV) + V • f . (100) 
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A.l - Physical/conditional observables - Hidden 
variables 



The fluid fields {Z} are, by assumption, prescribed 
smooth real functions of (r, t) S f2 X J. In particular, 
they can be either physical observables or conditional ob- 
servable, according to the definitions indicated below. 



Definition - Physical observable/conditional 
observable 



A physical observable is an arbitrary real-valued and 
uniquely-defined smooth real function of (r,t) 6 0x7. 
Hence, as a particular case (r,t) are observable too. 

A conditional observable is, instead, an arbitrary real- 
valued and uniquely-defined smooth real function of 
(r,i) eflxJ which depends also on non-observable vari- 
ables and is, as such, an uniquely-prescribed function of 
the latter ones. 

Therefore the functions Zi can be assumed respectively 
of the form [lj, [3] 



Zi = Zi(r,t) 



Zi(r,t,a), 



(101) 



(102) 



ct G V a C K fe (with k > 1) denoting a suitable set of 
hidden variables. In fluid dynamics these are intended 
as: 



Definition - Hidden variables 

A hidden variable is as an arbitrary real variable which 
is independent of (r, t) and is not an observable. 

A. 2 - Deterministic and stochastic fluid fields 



Hence, fluid fields of the type (|102[) are manifestly non- 
observables. However, if in the whole set flxl xV a , they 
are uniquely-prescribed functions of (r, t, a) then they are 
conditional observables. Hidden variables can be con- 
sidered in principle either deterministic or as stochastic 
variables, in the sense specified as follows. 



Definition - Stochastic variables 

Let {S, S, P) be a probability space; a measurable func- 
tion a :S — ► V a , where V a C R fc , is called stochastic (or 
random) variable. 

A stochastic variable a is called continuous if it is en- 
dowed with a stochastic model {g a , V a } , namely a real 



function g a (called as stochastic PDF) defined on the set 
V a and such that: 

1) g a is measurable, non-negative, and of the form 



9a{r,t, •); 



(103) 



2) if A C V a is an arbitrary Borelian subset of V a 
(written A £ B(V a )), the integral 



P a (A) = J dxg a (r,t,x) 



(104) 



exists and is the probability that a GA; in particular, 
since ex £ V a , g a admits the normalization 



J dxg a (r,t,x) = P a {V a ) = 1. 



(105) 



The set function P a : B(V a ) -> [0, 1] defined by (fT04f is 
a probability measure and is called distribution (or law) 
of a. Consequently, if a function f:V a — ► Vj C K m is 
measurable, / is a stochastic variable too. 

Finally define the stochastic- averaging operator 
OJsee also [HEl) as 



(f) a = (f(y, -))a = / dx 5-(r, t, x)/(y, x), (106) 

V a — * M, where 



for any P Q -integrable function /(y 
the vector y is some parameter. 



Definition - Homogeneous, stationary stochastic 
model 



The stochastic model {g a , V a } is denoted: 

a) homogeneous if g a is independent of r, namely 

b) stationary if g Q is independent oft, i.e., 

Sa = 0a(lV). (108) 

Definition - Deterministic variables 

Instead, if g a (r,t, ■) is a deterministic PDF, namely it 
is of the form 



? Q (r,t,x) =<5 (fc) (x-a ), 



(109) 



x — a ) denoting the fc-dimensional Dirac delta in 
the space V a , the hidden variables a are denoted as de- 
terministic. 

Let us now assume that, for a suitable stochastic model 
{9a, V a }, with g a non-deterministic, the stochastic vari- 
ables Zi = Zi(y,t,o) and /i(r, v,i, a) (where Zi(r,t,-) 
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and /i(r, v, t, •) are measurable functions) admit every- 
where in f2 x I and T x I the stochastic averages (Zi) a 
and (fx} a defined by (fTM)l . 

Hence, Zi = Zi(v,t,a), /i(r,v, t, a) and the mean- 
field force F(/i) [see Sections 2,3 and 4] admit also the 
stochastic decompositions 

Zi = (Z i } a + SZ i , (110) 
h = (fi) a + Sfx, (HI) 
F(A) = (F(/ 1 )) Q + ( 5F(/ 1 ). (112) 

In particular, unless g Q (r, t, •) is suitably smooth, it fol- 
lows that generally (Zi) , SZi and respectively , Sfi 
may belong to different functional classes with respect to 
the variables (r, t) . 



Then iVi(rj, v,i, a; 2T) can be defined as follows 



Nx(ri,v,t,a;Z) = 

f II e lfc (v)e ( ^ 

fe=l,2,3 



(114) 



M 
V B 



«*). (H5) 



e ifc ( v ) = e(Vfc(ri,t)-» fc -^:) (116) 

6(^-14^,*) + ^-), 

with M G N and 0(a:) the Heaviside theta function; here 
c G M and G N are defined so that there results 



v% e n e ^-( v ) 



f=l.JV fc=l,2,3 



M 3 = AT. 



(117) 
(118) 



A. 3 - Deterministic and stochastic INSE problems - 
Regular and turbulent flows 

Therefore, assuming, for definiteness, that all the fluid 
fields Z, the volume force f and the initial and boundary 
conditions, are either deterministic or stochastic variables 
and both belong to the same functional class, i.e., are 
suitably smooth w.r. to (r, t) and a, Eqs. (f88|) - [93|) define 
respectively a deterministic or stochastic initial-boundary 
value INSE problem. In both cases we shall assume that 
it admits a strong solution in f2 x / (or £1 x I x V a ). 

In the first case, which characterizes flows to be de- 
noted as regular, the fluid fields are by assumption physi- 
cal observables, i.e., uniquely-defined, smooth, real func- 
tions of (r,t) G f2 x I [with Q, the configuration space, 
and ft its closure, to be assumed subsets of the Euclidean 
space on R 3 and /, the time axis, denoting a subset of 
E]. 

In the second case, characterizing instead turbulent 
flows, the fluid fields are only conditional observables (see 
again Subsection A.l). In this case, besides (r,t), they 
may be assumed to depend also on a suitable stochastic 
variable a, (with a GV a and V a subset of R fc with k > 1). 
Hence they are stochastic variables too. 



APPENDIX B: DEFINITION OF Ni 



Let us define provide here an explicit definition of 
Ni(ri,v,t, a; Z) [required to specify also J^ req ^ in terms 
of Eq. (fT5)) ]. For definiteness, let us assume that the fluid 
velocity is bounded, i.e., that there exists Vg 6 M such 
that in Q x I for each component of the fluid velocity 
Vk(r,t, a) (with k = 1, 2, 3) there results 



Thanks to positions (|114[) - (|118p for an arbitrary N e N 
fulfilling Eq.((TTU), it follows 

i^SoTf zZ [ d 3 vN 1 (r l ,v,t,a;Z) = l. (119) 

A 1 AT *■* 



i=l, N 



Hence Eq. (jT9]) is satisfied identically. 

APPENDIX C: EVALUATION OF F (T) 

In /i(l) and /i(2) coincide with a local Gaussian [i.e., 
see Ea. (j55| ] there results by construction 

J rf 3 VV/ 1 (l)/ 1 (2) = \pv-U\f 2 {T U T 2 ,v;ZiV2Q) 

J d 3 W 2 /i(l)/i(2) = ^/ 2 (ri,r 2 ,v;Z), (121) 



where there results 

< P (1)"< P (2) 



(3 



< P (1)+< P (2) 
1 



(122) 



u= V(l)^(2) + V(2yi) = 

V(l) . V(2) 



v'f, (1) 1 v'i, (2) 

th,p\ l t h,p \ ' 

l I 

7V foT 



|^(r,t,a)| < -V B . 



(H3) 



Let us now evaluate the expression 
Jd 3 Vl[F(l)-F(2)]/ 1 (l)/ 1 (2). Introducing the 
notations 

F (l) = — f(l) + [V + v-V(l)]- ViV(l)+(124) 
Po 

+*A7?V(1), 
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Fi(1)= [V + v-V(l)] A(1)+ (125) 

+ ^) VlPl{1) \^JT) 2)' 



F (2) = — f(2) + [V-v- V(2)] • V 2 V(2)+(126) 
Po 



Fl(2)= [V-v -V(l)] 



(127) 



< P (2) 



V 2 pi(2) 



2pi(2) 
and letting for j = 1,2, 

m = 

it follows 



[V-v-V(2)] 2 3 
<^2) 2 



1 L> 



piU) 



(128) 



| d 3 vi[F(l)-F(2)]/ 1 (l)/ 1 (2) 
h = F (T) / 2 , 



Ff)+FP 



where 



2F 



(T) 



if(l)- 
Po 



+ [(/? + l)v-U-V(l)]-ViV(l) 

+i/V?V(l) - 
1 



-f(2) 



- [(/3 + 1) v - U - V(2)] • V 2 V(2)+i/V|V(2) 



2Ff ) = [(/? + D v-V(l)] A(1) , 



2pi(l) 



Hfc + 2 [/3v - U] ■ [v - V(l)] + [v - V(l)]^ _ 3 



[(/? + l)v-V(2)] 



A(2)- 



V 2 pi(2) 



<g(2) 
2pi(2) 

' f <4 + 2 [/3v - U] • [v - V(2)] + [v - V(2)] 2 3 ' 



< P (2) 



[1] E. Hopf, J. Ration. Mech. Anal. 1, 87 (1952). 



[2] G. Rosen, Phys. fluids 3,19 (1960). 
[3] S. E Edwards, J. Fluid Mech. 18, 239 (1964). 
[4] G. K. Batchelor, The Theory of Homogeneous Turbu- 
lence, (Cambridge University Press, Cambridge, Eng- 
land, 1960). 

[5] C. Dopazo, in Turbulent Reacting Flows (P. A. Libby and 
F.A. Williams Eds., London Academic Press), Chapter 
7, 375 (1994). 

[6] S.B. Pope, Turbulent flows, Cambridge University Press, 
p.463 (2000). 

[7] A. S. Monin, J. Appl. Math. Mech. 31, 1057 (1967). 
[8] T. S. Lundgren, Phys. Fluids 10, 969 (1967). 
[9] A.S. Monin and A.M. Yaglom, Statistical Fluid Mechan- 
ics, Vol. 1 and 2, MIT Press (1975). 
[10] S.B. Pope, Phys. Fluids 26, 3448 (1983). 
[11] M. Tessarotto and M. Ellero, RGD24 (Italy, July 10-16, 

2004), AIP Conf. Proceedings 762, 108 (2005). 
[12] M. Ellero and M. Tessarotto, Physica A 355, 233 (2005). 
[13] M. Tessarotto, M. Ellero and P. Nicolini, Proc. 26th 
RGD (International Symposium on Rarefied gas Dynam- 
ics, Kyoto, Japan, July 20-26, 2008); arXiv:0806.4530vl 
[physics.flu-dyn] (2008) . 
[14] M. Tessarotto, MHD 45, No. 2, 3 (2009). 
[15] I. Hosokawa, Phys. Rev.E 73, 067301 (2006). 
[16] G. Falkovich and V. Lebedev, Phys. Rev. Lett. 79 (21), 
4159 (1997). 

[17] Y. Li and C. Meneveau, Phys. Rev. Lett. 95, 184502 
(2005). 

[18] I. Hosokawa, Phys. Rev. E 78, 066312 (2008). 
[19] U. Frisch, Turbulence, (Cambridge University Press, 

Cambridge, U.K.), Chap. 3, 38 (1995). 
[20] M.J. Vishik and A.V. Fursikov, Mathematical problems 
of Statistical Hydrodynamics, (Kluwer, Dordrecht, NL, 
1988). 

[21] D. Ruelle, Chaotic Evolution and Strange Attractors, 

(Cambridge University Press, Cambridge, U.K., 1995). 
[22] M. Tessarotto and M. Ellero, Physica A 373, 142 (2007); 

arXiv:0806.4546vl . 
[23] M. Tessarotto, M. Ellero and P. Nicolini, Phys. Rev. A 

75, 012105,arXiv:quantum-ph/060691 (2007). 
[24] M. Tessarotto and M. Ellero, Proc. 25th RGD (Inter- 
national Symposium on Rarefied gas Dynamics, St. Pe- 
tersburg, Russia, July 21-28, 2006), Ed. M.S. Ivanov and 
(130) A.K. Rebrov (Novosibirsk Publ. House of the Siberian 

Branch of the Russian Academy of Sciences), p. 1001; 
arXiv:physics/0611113 (2007). 
[25] C. Cremaschini and M. Tessarotto, Proc. 26th RGD 
(International Symposium on Rarefied gas Dynamics, 
Kyoto, Japan, July 20-26, 2008); arXiv:0806.4546vl 
[physics.flu-dyn]; AIP Conf. Proc. 1084, 188 (2008). 
[26] M. Tessarotto, M. Ellero, and P. Nicolini, ibid.; 
arXiv:0806.4530v [physics.flu-dyn]; AIP Conf. Proc. 
1084, 33 (2008). 
[27] Marco Tessarotto, Claudio Cremaschini, Piero Nicol- 
ini, and Massimo Tessarotto, ibid.; arXiv:0806.4823vl 
[physics.flu-dyn]; AIP Conf. Proc. 1084, 182 (2008). 
[28] Marco Tessarotto and Massimo Tessarotto, ibid.; 
arXiv:0807.0274vl [physics.flu-dyn]; AIP Conf. Proc. 
1084, 483 (2008). 
[29] M. Tessarotto, M. Ellero, N. Asian, M. Mond, and P. 
Nicolini, ibid.; arXiv:0806.4817vl [physics.flu-dyn]; AIP 
Conf. Proc. 1084, 224 (2008). 
[30] M. Tessarotto, M. Ellero, D. Sarmah and P. Nicol- 



(129) 



21 



ini, ibid.; arXiv:0806.4931vl [physics.flu-dyn]; AIP Conf. 

Proc. 1084, 170 (2008). 
[31] E. Fonda, M. Tessarotto, P. Nicolini, and M. Ellero, 

ibid., E. Fonda, M. Tessarotto, P. Nicolini, and M. 

Ellero, ibid.; arXiv:0704.0339vl [physics.flu-dyn]; AIP 

Conf. Proc. 1084, 164 (2008). 
[32] Marco Tessarotto, Claudio Cremaschini, and Massimo 

Tessarotto, Physica A 388 (18), 3737-3744 (2009). 
[33] E.T. Jaynes, Phys. Rev. 106, 620 (1957). 



[34] A. Naert, R. Friedrich and J. Peinke, Phys. Review E 56 

(6), 6719 (1997). 
[35] R. Friedrich, J. Peinke, Physica D 102, 147 (1997). 
[36] St. Luck, J. Peinke and R. Friedrich, Phys. Rev. Lett. 

83, 5495 (1999). 
[37] Ch. Renner, J. Peinke, J. Fluid Mech. 433 383 (2001). 
[38] Ch. Renner, J. Peinke, R. Friedrich, O. Chanal and B. 

Chabaud, Phys. Rev. Lett. 89, 124502-1 (2002). 



